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ABSTRACT 

Two-photon decays of Higgs bosons are important channels for the search of 
these particles in the intermediate mass range at the pp coUiders LHC and SSC. 
Dynamical aspects of the Higgs coupling to two photons can also be studied by 
means of the 77 fusion of Higgs particles at high-energy e"'"e~ linear colliders. 
Extending earlier analyses which had been restricted to the Standard Model, 
we present in this note the QCD radiative corrections to the 77 decay widths 
of scalar and pseudoscalar Higgs particles in multi-doublet extensions of the 
Higgs sector, as realized for instance in supersymmetric theories. 
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The exploration of the electroweak symmetry breaking mechanism is one of the most im- 
portant tasks of particle physics. Within the Standard Model, the solution of this problem 
is associated with the existence of fundamental scalar Higgs particles. If the Higgs particle 
in the Standard Model [SAi] is light, with a mass in the intermediate range below ~ 200 
GeV, the theory can be extrapolated perturbatively up to the GUT scale. However, the 
hierarchy problem suggests the supersymmetric extension of the model in this case, sta- 
bilizing the mass of light Higgs bosons in the background of high-energy GUT scales. In 
contrast to the SAi, supersymmetric theories incorporate a spectrum of Higgs particles 
$, in the minimal version [A^iSiSAI] light and heavy scalar [CP-even] particle^ h, H, a 
pseudoscalar [CP-odd] particle A, and a pair of charged Higgs particles. 

The precise prediction of the 77 decay widths of Higgs particles [Q is important for two 
reasons. First, the 77 decay mode plays a crucial role for the search of these particles 
in the lower part of the intermediate mass range at the pp colliders LHC and SSC |@]. 
Second, the 77 widths can be measured directly by means of 77 fusion at high-energy 
e~^e~ linear colliders f^, ^. Since the photons couple to Higgs bosons via heavy particle 
loops, the 77 widths are sensitive to particle masses, standard and also supersymmetric, 
well above the Higgs masses themselves. However, these effects are small if the particle 
masses are not generated through the Higgs mechanism. To exploit this method, it is 
therefore mandatory to control properly the radiative corrections to the $77 couplings 
mediated by the standard particles. 



In this note, we shall present the QCD radiative corrections of the quark-loop contribu- 
tions to the 77 widths of the Higgs particles in the [minimal] supersymmetric extension 
of the Standard Model, Fig. 1. This demands the extension of the calculation for the SAi 
Higgs decay [^, in two ways: (i) Since the decay amplitudes in part of the SlASy pa- 
rameter space are dominated by 6-quark loops, the QCD corrections are to be analyzed 
for scalar h, H masses above the fermion-antifermion threshold^; and (ii) the QCD cor- 
rections are to be determined for pseudoscalar Higgs bosons A. 



The 77 couplings to Higgs bosons are mediated by charged heavy particle loops \W , 
fermion, chargino, sfermion and charged Higgs boson loops in the scalar case /i, if, and 
fermion and chargino loops in the pseudoscalar case A\. Denoting the quark amplitudes 
by Aq etc., the 77 decay rates are given by 0, ^ 



V{H - 77) 

V{A ^ 77) 
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^The scalar particles will generically be denoted by H. 

^This problem has first been solved by one of the authors The analysis would agree with the 
results of a recent preprint Q if the large logarithms in the QCD corrections were not mapped into the 
running quark mass [see Appendix A] . 
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where the quark and W amphtudes read at lowest order 
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The scahng variable is defined as r = M^/Amj with rrii denoting the loop mass, and 
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for r < 1 
for r > 1 
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The coefficients gf denote the couplings of the Higgs bosons [normalized to the SM. Higgs 
couplings] to top, bottom quarks and W bosons, recollected for the sake of convenience 
in Table 1. 



Table 1: Coefficients gf 
of the Higgs couplings to 
quarks and W bosons in the 
M.SSM.. a, (3 are mixing 
angles, tg/5 = V2/V1 be- 
ing the ratio of the Higgs 
vacuum expectation values. 
Numerical values of the co- 
efficients are presented in 
Ref. ITUi, for instance. 
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The cross section for the 77 fusion of Higgs bosons is found by folding the parton cross 
section with the 77 luminosity, see e.g. 0. The parton cross section is determined by 
the 77 width so that ((1(77 — » $)) oc r($ 77)/M|. The photon beams are generated 
either automatically as initial-state Weizsacker- Williams and beamstrahl photons in 
collisions, or they may be generated, in a dedicated effort, by Compton back-scattering 
of laser light p[ . 



The QCD corrections to the quark amplitudes can be parametrized as 
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LO 
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The coefficient C depends on r = M^/AmQlfi^) where the running quark mass mQ{ii'^) 
is defined at the renormalization point /i which is taken to be /i = M$/2 in our analysis; 
this value is related to the pole mass mQlmq) = mq in the on-shell renormalization 
scheme by 



(6) 
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The lowest order amplitude Aq-^ is to be evaluated for the same mass value mQ([M$/2]^). 
The choice /i = Mq,/2 of the renormalization point avoids large logarithms logM|/mg 
in the final results for Higgs masses much larger than the quark mass [for details see 
Appendix A], is taken at /i for A = 200 MeV. 



We have evaluated the diagrams of the type shown in Fig. 1 plus the corresponding counter 
terms for the running quark mass at /i = M$/2. The 't Hooft-Veltman 75 prescription 



has been adopted for the dimensional regularization of the amplitudes. The 5-dimensional 
Feynman parameter integrals have been reduced analytically to 1-dimensional integrals 
which have been calculated numerically. 



The amplitudes C-j^ for scalar loops and Ca for pseudoscalar loops are shown in Fig. 2a/b 
as functions of r. The coefficients are real below the quark threshold r < 1, and complex 
above. Very close to the threshold, within a margin of a few GeV, the present perturbative 
analysis^ cannot be applied anymore. [It may account to some extent for resonance effects 
in a global way.] Since QQ pairs cannot form 0"'""'" states at the threshold, XmC-n vanishes 
there; TZeC-n develops a maximum very close to the threshold. By contrast, since QQ 
pairs do form O""*" states, the imaginary part XitiCa develops a step which is built-up 
by the Coulombic gluon exchange [familiar from the Sommerfeld singularity of the QCD 
correction to QQ production in e^e~ annihilation]; TZcCa is singular at the threshold. For 
large r, both coefficients approach a common numerical value, as expected from chiral 
invariance in this limit. In the opposite limit, the QCD corrections can be evaluated 
analytically, 

mg > M$ : Cu ^ -1 and Ca ^ 

These results can easily be traced back to the form of the 77 anomaly in the trace of the 
energy-momentum tensor |jl3[ and to the non-renormalization of the axial-vector anomaly 



1^, as demonstrated in Appendix B. 



In Fig. 3a-d the QCD corrected 77 widths for H, A Higgs bosons are displayed in the 
minimal supersymmetric extension of the Standard Model [taking into account only quark 
and W boson loops] for two values tgP = 2.5 and tg/9 = 20. While in the first case top 
loops give a significant contribution, bottom loops are the dominant component for large 
tgjS. The overall QCD corrections are shown in the lower part of the figures. The correc- 
tions to the widths are small, ~ 0{as/Tr) everywhere. [Artificially large 6 values occur 
only for specific large Higgs masses when the lowest order amplitudes vanish accidentially 
as a consequence of the destructive interference between W and quark-loop amplitudes, 
see also [|, |.] 

In conclusion. We have calculated the QCD corrections to the decays h,H,A^ 77 of the 
neutral scalar and pseudoscalar Higgs bosons in the minimal supersymmetric extension 
of the Standard Model, that is taken for illustration. These corrections are well under 

^By choosing the renormalization point ^ — M$/2 the perturbative threshold Eth — 2mQ{mQ) 
coincides with the on- mass shell value proper. A shift between M$/2 and M$, for instance, affects the 
widths very little away from the threshold. 
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control across the physically interesting mass ranges, if the running of the quark masses 
is properly taken into account. 



APPENDIX A 

The QCD corrected quark amplitude for $ ^ 77 may be written in the general form 



AQ = 4^(mQ)<'l + 



ciimq) + C2(mQ)log — 



71 



where rriQ = mQ{^^) is the quark mass defined at the renormalization point /i. The scale 
in Us may in principle be chosen different from /i. For large mg, ci approaches —1 for 
/i, if, and for A, while C2 vanishes ~ l/roQ, and no large logarithms appear in this limit. 
For niQ —>■ 0, however, C2 approaches a finite non-zero value while ci develops a large 
logarithm the coefficient of which is given by C2 , 

C2 ^ 2 

Ci — C2 log (M|/4mQ) + const 

By choosing fi = M<j,/2, all large logarithms are eliminated from the coefficient of as and 
mapped into the effective quark mass of the lowest-order amplitude. This is reminiscent 
of the corresponding procedure for Higgs decays to fermions pairs [|16|. Had we chosen 
/i = niQ instead 0, we would be left with unnaturally large corrections not taking ad- 
vantage of renormalization group improvements. 

A technical remark ought to be added on a subtle problem related to the 't Hooft-Veltman 
implementation of 75 in the dimensional regularization scheme which reproduces the axial- 
vector anomaly to lowest order []T2| automatically. The multiplicative renormalization 



factor of the scalar QQ current is well-known to be given by Zt-^qq = 1 — Z2 + Smq/mQ 
where Z2 is the wave-function renormalization factor and dniq the additive mass shift. 
To insure the chiral symmetry relation T^{p\p) '-y^T{p\p) in the limit itlq for 
the fermionic matrix element of the pseudoscalar and scalar currents, the renormalization 



factor of the pseudoscalar current has to be chosen [jT7 



Zaqq = Z-HQQ + 8as/(37r) 



The additional term supplementing the naive expectation is caused by spurious anomalous 
contributions that must be substracted by hand [jTSj . 
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APPENDIX B 



For large quark masses compared to the Higgs mass mg 3> M$, low energy theorems can 
be exploited to calculate the corrections C-j-i^a in this limit. 



The scalar coupling 7-^77 can be derived |19| from the requirement that the matrix element 



of the trace of the QCD corrected energy-momentum tensor [jT3 



vanish in the low energy limit; 5 = 2as/7i and = 2a^/7r{l+as/7r+- ■ ■) is the QED/QCD 
/3-function. Since the Higgs bosons are coupled to the mass operator (mg/f) QQ, the 
QCD corrected effective Lagrangian can be determined immediately, 

>Ceff(H77) = ^ {V2G,f\l (1 - ^) HF.^F^'^. 



From the non-renormalization of the anomaly of the axial-vector current [|1^, [T^ follows 
the non-renormalization of the A77 coupling, 

d^A^ = 2mQQtj5Q + -tF^uF^" 

4:71 

SO that the effective Lgrangian 

£efr(A77) = ^ (V2GFy^' AF^^^ 
is valid to all orders of as for the two-photon irreducible part of the diagrams. 



Note added. As a corollary to this result we observe that the irreducible part of the 
Agg coupling will not be renormalized either, and only the reducible diagrams need be 
calculated for the QCD corrections to the production process gg ^ A in the range 
iTiA < 2mQ. Only the virtual corrections are different from the scalar case, C 

+ 6 + 1(33 - 2 Aj.) log /iVm^ in the notation of eq.(9) in ReL'^. [A. D., Madi- 
son SSC Workshop '93; see also R. Kauffman and W. Schaffer, BNL preprint).] 
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FIGURE CAPTIONS 

Fig.l. Generic diagram of the QCD radiative corrections to the Higgs coupUng to two 
photons. 

Fig.2. (a) Real and imaginary part of the radiative corrections to the quark amphtude 
of the scalar 7-^77 coupling, normalized to the lowest-order amplitude; (b) the same 
for the pseudoscalar ^77 coupling. 

Fig. 3. (a) Two-photon widths of the AiSSAi Higgs bosons h, H, A for tgf3 — 2.5, and 
(b) size of the QCD radiative corrections; (c), (d) the same for tg(3 — 20. 
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